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It is shown that photon-photino inter-conversions (if exist) 
may result in a detectable CMBR spectrum distortion which 
amplitude depends on photino properties, such as its mass. 
An upper bound on the distortion parameter determined from 
the recent COBE-FIRAS data, hence, sets a lower bound on 
the photino mass, m 7 > 300 eV. 
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The elegant and generally accepted theory of super- 
symmetry in particle physics (for astrophysical implica- 
tions, see, e.g., which unifies bosons and fermions 
remains neither confirmed nor rejected experimentally 
so far. This theory predicts the existence of new par- 
ticles, the super-symmetric partners for all known par- 
ticles, neither of which has been detected yet. In this 
Letter, we show that, under particular conditions, the 
photino 7 (the ^-spin fermionic super-symmetric partner 
of the photon 7), if exists, may result in a measurable dis- 
tortion of the blackbody Cosmic Microwave Background 
Radiation (CMBR) energy spectrum. The magnitude of 
the deviation from the Planckian spectrum is directly re- 
lated to the value of the photino mass, to 7 ^ 0. We then 
use recent results j| from the COsmic Background Ex- 
plorer (COBE) mission to put some constraints on the 
value of m~. Recent evidence for the time variability of 
physical constants (e.g., the fine structure constant) ||] 
makes the measurement of the 'photino induced' CMBR 
spectrum distortion to be of great importance. Indeed, 
the CMBR provides information on physical laws existed 
in the early Universe (at redshifts z ~ 10 3 — 10 8 ) which 
cannot be obtained by means of other astronomical ob- 
servations Q. Thus, such an effect, if detected, would 
have profound implications for our understanding of fun- 
damental physics. 

Photons during the radiation dominated era (10 3 < 
z < 10 9 ) are in thermal and ionization equilibrium with 
matter. However, the process of photon energy equili- 
bration throughout the spectrum is quite slow, because 
it proceeds via the (inverse) Compton scattering with 
electrons Q , which number density is small with respect 
to the number density of photons, n e /n 7 ~ 10 -9 . In 
fact, the rate of Comptonization of photons becomes 
comparable to the Universe expansion rate at redshift 
z c ~ 2 x 10 4 r2 ~ 10 5 (where f^o = po/p cr it is the 
ratio of present density of matter in the Universe to the 
critical density). Thus, the characteristic time, T eq , re- 
quired to establish the (spectral) equilibrium in the pho- 
ton gas diverges as z approaches z c . 



On the other hand, in the super-symmetry theory, 
the photon and photino arc different "spin states" (with 
s = 1 and s = i, respectively) of a particle and have 
different interchange properties. The photon is known 
to be a massless particle, and the photino is believed 
to be massive. Thus, if a corresponding mixing angle 
8 of these two 'states' is not identically zero, one may 
expect the 7^7 oscillations to occur, with the prob- 
ability P°£^(i) = § sin 2 29 [1 - cos ((J% - E^)trH)] , 
where the energy difference is — E 7 ~ m 7 c 2 . Such 
a simple picture, however, is not correct, because the os- 
cillations of a single particle are forbidden due to angular 
momentum and R-parity conservation. Thus, the 7 «-> 7 
inter-conversions are allowed to occur only in reactions 
(collisions) with other particles. (We postpone the dis- 
cussion of probable reaction channels intill the end of 
the Letter.) Consequently, the average rate z^ 77 of 7 <-> 7 
inter-conversions is, thus, the average collision frequency, 
(<7^v)n, times the probability of a 7 <-> 7 oscillation dur- 
ing collision time, P° s ^(At cM ), (note, P^- = P?^, 
by symmetry). It should be emphasized that 7 <-> 7 
inter-conversions are incoherent in this regime, i.e., the 
quantum coherence of photon-photino oscillations is com- 
pletely lost at times greater than v~- due to classical, in- 
trinsically chaotic motion of interacting particles. Photi- 
nos are, hence, in (at least, marginal) equilibrium with 
photons. The equilibrium fraction of 7's in a system is 
proportional to the forward-to-reverse transition rate ra- 
tio and, hence, is independent of the mixing angle. The 
rate with which equilibrium is established is, however, 
affected by 0. We now introduce the total characteristic 
time of a 7 <-> 7 transition as follows, 
The numbers of photons and photinos in the system, 
thus, fluctuate on the time-scale r 77 , while their total 
number is preserved. 

Let us consider two limiting cases. First, the spec- 
tral equilibration time is much smaller than the typi- 
cal 7 <-> 7 fluctuation time, r eq -C t 77 . This is a gen- 
eral case, because super-symmetric particles are believed 
to be weakly interacting (otherwise they would be de- 
tected). In this case, the number of photons in the sys- 
tem fluctuates adiabatically slowly and CMBR photons 
have enough time to re-arrange throughout the spectrum. 
Hence, the Planckian (or Bose- Einstein H) spectrum is 
maintained. Second, if the spectral thermalization pro- 
cess is (externally) inhibited, e.g., by the Universe expan- 
sion, the opposite case may realize. Then 

T 77 <C T eq , (l a ) 
^77 texpan ^ (3-P) , (^^) 
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(where H is the instantaneous Hubble constant) and 
7 <-> 7 transitions are faster than the photon transit time 
though the spectrum. Thus, different, 'mean-field' equi- 
librium, which is, in general, a weighted sum over all 
realizations, is created; rather than the 'instantaneous' 
Planckian equilibrium. [Note, inequality ( |Tb| ) ensures 
that 7 <-> 7 inter-conversions are fast enough, compared 
to the expansion rate]. The partition function and oc- 
cupation numbers for a gas of particles with fluctuating 
quantum statistics were calculated elsewhere ||. 

Since the CMBR spectrum is very close to Planckian 
with temperature To ~ 2.7 K, we look for a small devia- 
tion only. Then, two free parameters, pf and pb (the bo- 
son and fermion probabilities, see Ref. Q for details), be- 
come pt = 1—pf and pf <C 1. Thus, we may expand the 
occupation number n in terms of pf. From Eqs. (11,12) 
of Ref. lU?]], we straightforwardly obtain the effective 
bosonic and fermionic contributions (rib + rl / = n): 



n 7 — rib 




yfe x (e x - 1) ' 
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where x — e/kT is the dimensionless energy. Thus, the 
deviation of the photon distribution from the Planckian 
one is 



Sn 7 = n 7 — npi = y/pj 



( e x _ 1)3 



(3) 



The parameter pf is related to the photino mass, m 7 , 
as follows. Let us notice that condition (|l^) states 
that photinos are in equilibrium with photons: c>tn 7 = 
Pb^fR-y — -P/^/,n 7 = 0. The probability to change a 
'state' from bose to fermi, Pb->f, is complimentary to that 
to stay in a bose 'state', pb- Hence Pb^f = 1 — Pb = Pf, 
and we obtain: 



Pf — — = — 
Pb n 7 



(4) 



The number densities of photons and photinos are cal- 
culated integrating Eqs. (||) over energies. Photons are 
relativistic (e = pc) and the degeneracy g 7 = 2 (due to 
two polarizations); photinos are assumed non-relativistic 
(m 7 c 2 ^> kT c ) and <? 7 = (2s + 1) = 2. Finally, we obtain: 
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C(3) V 8 V kT, 
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-m=fC /kT c 



(5) 



1.2 and T c is 
2.7 x 10 5 K, 



where £ is the Riemann ^-function, £(3) 
the temperature at redshift z c , T c ~ T$z, 
to satisfy condition (|Ta|). 

The CMBR spectrum has been measured with great 
accuracy by the FIRAS (Far-InfraRed Absolute Spec- 
trophotometer) on board of the COBE satellite. For our 



analysis, we use the CMBR residuals (after subtraction 
of the Planck blackbody and galactic emission spectra) 
of both monopole and dipole components of Ref. [0 (re- 
ferred to as the '96-monopole' and '96-dipole' data sets, 
respectively) and, for comparison, some older data of Ref. 
g ('94-monopole' data set) and Ref. @ ('94-dipole' data 
set). Note, if a dipole component is associated with the 
Earth motion only (with respect to the CMBR), both 
monopole and dipole spectra should be identical. In gen- 
eral, a dipole spectrum is less sensitive to systematic er- 
rors (e.g., absolute calibration), but statistical errors may 
be higher, because a dipole component is of much smaller 
amplitude than monopole. 

Since deviations from the Planck spectrum are small, 
a linear fit can be performed: 



J2 a 

models(a) 



dSgjl 

da 



(6) 



Here B v is the Planck blackbody spectrum, S a repre- 
sents a spectral model, and a is a distortion parameter 
which quantify the deviation from B v (Tq). It is impor- 
tant to have the first term since dB v jdT correlates with 
other models dS a /da. We consider three spectral mod- 
els. First, the 'photino induced' distortion (referred to as 
the Jpf- model): 



dS 



dVPJ 



= 2/icV 




(7) 



which readily follows from Eq. (^5j), where xq = hcu/kTQ 
and v is measured in cm" 1 . The distortion parameter is 
y/pj- Second, the Bose-Einstein quasi-equilibrium pho- 
ton distribution H with dimensionless chemical potential 
fi (ji- model). Third, the Comptonized spectrum charac- 
terized by the Kompaneets [|lfj parameter y (y- model). 
(For details regarding the last two models, see, e.g., Ref. 
|g|.) Note that all three distortion parameters, y/pj, H, 
and y, must be positive by physical meaning. The nor- 
malized spectrum distortions and the effect of a tem- 
perature shift dB v /8T are shown in Fig. |]. Note very 
strong similarity of the y/pj- and /i-models, which pre- 
cludes their simultaneous fit. However, preference to one 
of the models may be given, due to their opposite contri- 
butions to the CMBR spectrum. 

Fits (||) have been performed with the help of Math- 
ematica 3.01 standard package for Linear Regres- 
sion analysis, which uses the weighted least square 
E l/ (A 1 ,/ov) 2 /X^( 1 / c 'v) 2 ] method, where A„ is the dif- 
ference of the model and the data at frequency u, and 
a~ v is the uncertainty. Test fits of the y-model and the fi- 
model separately to the both 96- and 94-monopole spec- 
tra yield values of y and /i consistent with those of Refs. 
pLpL I n ^ ne dipole fits, y and \x tend to negative (unphys- 
ical) values, so do their simultaneous fits, i.e., the (jj, + y)- 
model, for both monopole and dipole components. All 
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FIG. 1. The normalized spectrum dis- 

tortions (dSc / da) / \(dS a /da)\ max for the y/pj-model (bold 
solid line), /x- model (dashed line), y- model (dotted line), and 
the effect of a temperature shift (dB v / dT) / \(8B V / dT)\ max 
(solid line) vs. frequency v. 

TABLE I. Results of fits for Jp]- and {y/pj + j/)-models. 



Data set 


y/py-model 


(yjp] + y)-model 


^pjx 10 5 




y x 10 6 


96-monopole 


1±4 


3±8 


2±6 


96-dipole 


13±2 


26±4 


11±3 


94-monopole 


10±9 


50±14 


35±11 


94-dipole 


7±5 


3±9 


-4±6 



the fits, which include the both fi- and ^/py-models, yield 
large statistical errors resulting from very strong correla- 
tion of these models, so they have been rejected. Thus, 
there are only two statistically reliable models, namely 
the ^/py-model and the (y/pj + ?/)-model. The results 
with la uncertainties (95% confidence level) are pre- 
sented in Table H. There is an obvious inconsistency be- 
tween the 96-monopole and 96-dipole fits for both mod- 
els. Since the results for the 96-dipole and 94-monopole 
spectra calculated in the ^/py-model are also inconsis- 
tent with those in the {^fpj + y)-model (even within 3<r 
uncertainty), we consider them as statistically unreliable. 
What causes such anomalously large correlations of resid- 
uals (especially in the 96-dipole data set) is not known. 
This issue requires further investigation. (One of possi- 
ble 'candidates' is residual galactic contamination, which 
is strongly coupled to the y-model.) Other fits, i.e., for 
the 96-monopole and 94-dipole data, with both models 
yield results which are consistent with each other and in- 
clude the 'null' result (^/pj = y = 0) within ~ 1.4er un- 
certainty. Weak trend toward positive values may, how- 
ever, be discerned. Using the values with smaller statisti- 
cal uncertainties, we set the following conservative upper 
limit on the parameter ~Jp]'. 



VpI<^~ 4 . (8) 

From Eq. (|5|), a lower bound on the photino mass is 



which is much lower Q than generally accepted: m 7 ~ 
0.1 - 10 GeV/c 2 . Note, however, that the 'null' result 
cannot yield absolutely reliable lower bound on to 7 , be- 
cause such a result can also be due to breaking of condi- 
tion (lb), which accounts for 7 <-* 7 kinetics. 

Now we check the consistency of our result with con- 
ditions (1). In fact, we should check ( |Tb| ) only because 
( |la| ) is satisfied at z ~ z c by definition. Since binary col- 
lisions are most frequent, the leading processes are: (i) 
the binary photino collisions, 77 ^ 77; (ii) the R-parity 
'exchange', 7X 72;, where x is any particle (possibly an 
electron, e, as the most abundant) and 'tilde' denotes its 
super-symmetric partner; and (iii) the decay-type pro- 
cesses, 7 <=^ £7, where x may be an axino, a (the super- 
symmetric partner of the axion) , or another hypothetical 
particle. We estimate cross sections of these processes 
and show that only the last reaction cannot be com- 
pletely ruled out. Of course, the reaction channels that 
could influence the 7^7 transitions are not limited to 
those written above and may include, for instance, 'chain 
reactions', 7x1 — > 7x2, 72:3 — > 7x4, etc. Since complete 
study of 7 «-> 7 reactions is beyond the scope of this Let- 
ter, we restrict ourselves considering the transition pro- 
cesses mentioned above. Since the results below are il- 
lustrative, we assume, for simplicity, that sin 2 29 ~ 1 and 
(Ej — E 7 ) ~ to 7 c 2 3> hAt~^ u , so that the fast oscillating 
cosine term vanishes upon time and ensemble averaging. 
However, the mixing angle and oscillation time contribu- 
tions may be trivially recovered when needed. 

First, let us consider 77 77 process. In equilibrium, 
according to ([it]), (erv^n 2 = (cru)~~n?, where a is the 
effective cross-section of a process. The characteristic 
times are t_7 - ~ (av) n-» and riT „ ~ (cru)--n^,, and 

©] 



. ~~- (av) and t- *„ ~ (crv)-~n~, 

7 — > 7 \ I 77 J 7 — > 7 \ / 77 / 3 

3H may be replaced [due to condition (|Ta|)] with the 
Compton scattering inverse rate r" 1 ~ (a^vl h )n e , where 
cr c is the Compton (Thompson) cross-section and t>* /l is 
the electron thermal velocity. Then, from the definition 
of t 77 and ( |lb| ) and using pj ~ n 7 /n 7 , we express the 
effective cross-sections as follows: 



^77 > l^e _ 1q7 

(<r c vf) ~ pj n 7 



(a c vf) ~ n 7 



Hence, the 77 cross-section is too large compared to that 
expected pd| ] from particle physics experiments. [Such 
a cross-section corresponds to the photino mass (using 
@) m 7 < 0.5 eV/c 2 , which contradicts to @.] Thus, 
such a process is (probably) ruled out. Second, the same 
is true for jx "fx process. Indeed, the cross-sections 
(estimated for electrons, x = e) are 



(av) - 
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(<™>2g > ne 



„th\ 



(11) 



which are also too large. Third, the situation is better 



m 7 > 300 eV/c 



(9) 



for 7 <^ 07 process. Denoting 
rate, we obtain: 



decay 



to be the 7-decay 
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-pdecay ^ 
7 ~ 



3H C 
Pf 
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" 2 sec" 1 



ns, 



(12) 



where _ff c = H(z c ). Since the axino is a possible non- 
barionic dark matter candidate, the ratio n e /ng may be 
quite small, too. Estimated from the 7-decay rate [ p^[ , 
the photino mass is > 3 MeV/c 2 (for the Peccei- 
Quinn symmetry breaking scale Vpq — 10 9 GeV). Since 



yC 2 3> kT c , fine-tuning 



ma is required; otherwise 



in 

the inverse decay 0,7 — > 7 is strongly suppressed. Assum- 
ing m^, ~ rria ~ 3 MeV/c 2 , we may estimate the cosmo- 
logical abundance of axinos (for the fltot — 1 Universe) 
to be na/n e ~ ^Q 1 'nib ar i on /ma < 10 5 , which yields a 
reasonable 07 cross-section. 

It is important to note that the restrictions imposed 
by condition (|l^) may, however, be significantly relaxed 
by admitting (i) multi-channel 7^7 inter-conversions 
(e.g., 7^7 and 7 — > 7 conversions may proceed through 
different processes), (ii) super-symmetry theories beyond 
the standard model, or (iii) theories with R-parity viola- 
tion. Detailed study of these issues is, obviously, beyond 
the scope of this Letter. 

To conclude, we have shown that, under appropri- 
ate conditions, the 7 <-> 7 inter-conversions may result in 
a detectable distortion of the CMBR blackbody energy 
spectrum. Since the 'bare' 7^7 process is forbidden by 
angular momentum and R-parity conservation, kinetics 
of super-symmetric particle-matter interactions plays an 
important role. Thus, the amplitude of the CMBR spec- 
trum distortion is related to photino properties, e.g., its 
mass, rrij. Recent data from the COBE- mission are used 
to evaluate the distortion parameter y/pj and photino 
mass. 
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